Medium Modifications of the p Meson in Nuclear Photoproduction 
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We extend our recent study of dilepton invariant-mass spectra from the decays of p mesons 
produced by photon reactions off nuclei. We specifically focus on experimental spectra as recently 
measured by the CLAS Collaboration at the Thomas Jefferson National Accelerator Facility using 
carbon and iron nuclei. Building on our earlier work, we broaden our description to a larger set of 
observables to identify sensitivities to the medium effects predicted by microscopic calculations of 
the p spectral function. We compute mass spectra for several target nuclei and study the spectral 
shape as a function of the 3-momentum of the outgoing lepton pair. We also compute the so-called 
nuclear transparency ratio which provides an alternative means (and thus consistency check) of 
estimating the p width in the cold nuclear medium. 
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I. INTRODUCTION 



The study of in-medium properties of vector mesons 
has been pursued vigorously in recent years. The restora- 
tion of the spontaneously broken chiral symmetry at high 
temperatures and/or density predicts that hadron prop- 
erties change when approaching pertinent phase changes: 
the spectral functions of chiral partners are believed to 
degenerate [J-Q- Precursor effects of chiral restoration 
may already occur in more dilute hadronic matter vifhere 
medium modifications can be computed in more con- 
trolled calculations using many-body techniques. How 
far such medium effects can be related to fundamental 
properties of the QCD phase transition(s) remains an 
open question. On the experimental side, electromag- 
netic probes (real and virtual photons) are very suitable 
because no significant distortions occur owing to strong 
initial and/or final-state interactions with the nuclear 
medium. Invariant-mass spectra of dileptons in heavy- 
ion collisions (HICs) at the Super Proton Synchrotron 
(SPS) [B-Q clearly revealed the presence of medium mod- 
ifications of the p-meson spectral function in the hot and 
dense medium. The most recent NA60 dimuon data Q 
enable average medium effect on the p meson to be quan- 
tified as an approximately 3-fold broadening relative to 
its free width with little (if any) mass shift 0, [l^ . This 
average, however, results from a time evolution of a hot 
and dense fireball encompassing hadron densities which 
vary by approximately a factor of 5 from the critical 
temperature (Tc ~ 175 MeV) down to thermal freezeout 
(Tfo ~ 120-130 MeV). Thus, the predictive power under 
such conditions requires not only a good knowledge of the 
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time evolution of temperature and baryon density, but 
also a well-constrained model for the in-medium physics 
figuring into the vector-meson spectral functions. More 
elementary reactions using a single-particle projectile di- 
rected on nuclear targets offer a valuable simplification of 
both aspects: (a) the nucleus provides a static medium 
(at zero temperature) and (b) the medium consists of nu- 
cleons only (rather than including a tower of mesonic and 
baryonic excitations). A further benefit is that medium 
effects induced by nucleons (baryons) were indeed identi- 
fied as the more relevant one relative to mesons 
Although these features should, in principle, allow for a 
better theoretical control in the evaluation of medium ef- 
fects, there is, of course, a price to pay: (i) the medium 
density is restricted to that of normal nuclear matter, 
Qq = 0.16 fm""^, with significant density gradients ow- 
ing to surface effects, and (ii) the production mechanism 
is more involved (generally depending on the projectile) 
compared to an approximately thermal medium in ultra- 
rclativistic HICs. Nevertheless, if the latter point can 
be addressed satisfactorily (e.g., by gauging the produc- 
tion process on proton targets), valuable information on 
medium effects is expected to emerge from the investiga- 
tion of nuclear dilepton production experiments. 

Dilepton spectra in elementary reactions were mea- 
sured using both proton- and photon-induced production 
off nuclei. Using a 12 GeV proton beam at KEK, the 
E325 Collaboration reported a significant reduction in 
the p mass [l^l- However, with a 1-3.5 GeV incident- 
energy photon beam at JLab, the CLAS Collabora- 
tion [1^ [l^ , using an absolutely normalized background 
subtraction procedure, found not a significant mass shift 
but rather a moderate broadening of the dilepton excess 
spectrum associated with p decays (cf. also Ref. [Ol for 
two-pion production experiments on light nuclei). 

As indicated above, the starting point for a theoret- 
ical description of these reactions is a realistic model 
for the elementary production process, N e~^e~N 
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as studied by several authors [18l - l2l| . In our previous 
study [2 21. we combined the meson-exchange model of 
Ref. [2l[ with an in-medium p propagator [ifl] which de- 
scribes low-mass dilepton and photon spectra in HICs 
at the CERN-SPS 0. The production amplitude was 
augmented by baryon resonance channels to ensure con- 
sistency with the medium effects in the p propagator. 
Using an average local-density approximation for the de- 
cay points of the p, this model was found to describe 
the dilepton invariant-mass spectra off iron targets, as 
measured by the CLAS Collaboration, fairly well. 

In this article, we improve and extend our previ- 
ous work [i^l in several respects. First, we refine our 
schematic model for estimating the density probed by 
the p meson by folding over a realistic density distribu- 
tion (rather than using an average density), thereby ac- 
counting for momentum and density dependence of the 
decay rate. The improved description is compared to ex- 
isting carbon and iron data, supplemented by predictions 
for several heavier targets to illustrate the sensitivity to 
medium effects when going up to uranium. We further- 
more investigate the role of 3-momentum cuts on the 
outgoing lepton pair as an additional means to enhance 
the observable medium modifications. Finally, we com- 
pute the so-called nuclear transparency ratio as an inde- 
pendent (and thus complementary) observable to obtain 
quantitative information on the in-medium width of the 
p meson, by measuring how the absolute magnitude of 
the cross section depends on the nuclear mass number. 

This paper is organized along the aforementioned lines. 
In Sec.|ll]we briefly recall the main ingredients of our ap- 
proach, the improvements and the resulting density pro- 
files probed by the p decays. In Sec. IIIII we update the 
comparison to existing data (Sec. lIII A| and carry out the 
sensitivity studies with respect to the target (Sec. IIIIB[) 
and 3-momentum (Sec. IIII C|) dependence of the spec- 
tra, as well as the transparency ratio of the total yields 
(SeciniDl). We conclude in Scction|lVl 



II. p PRODUCTION AND DECAY IN NUCLEI 

We first recapitulate the main components of our 
model as constructed in Ref. [l^l- The elementary pho- 
toproduction amplitude, for the process 7p -> e+e^p, 
is based on the meson-exchange model developed by two 
of us (2]| . which accounts for cr, /2, it, and Pomeron 
t-channel exchange as well as nucleon s- and u-channel 
pole diagrams. This model gives a very good descrip- 
tion of the total p production cross section for incom- 
ing photon energies > 2.5 GcV. Since the input 
spectrum used by the CLAS Collaboration at JLab con- 
tains photon energies down to ~1 GeV, we extended the 
model by including s-channel baryon-resonance excita- 
tions. By implementing the same set of resonances (with 
the same coupling constants and form factors) that were 
used in the calculation of in-medium p spectral func- 
tion [13, , a good description of the low-energy part of 



the p-production cross section was found. The descrip- 
tion of the angular distributions (which are forward dom- 
inated for i-channel exchanges) also improves because the 
resonance decays generate significant strength at large 
scattering angles. This provides a realistic starting point 
for the production mechanism on nuclear targets and es- 
tablishes consistency between the production amplitude 
and the medium effects induced by the subsequent prop- 
agation of the p meson through the nucleus. In addition 
to direct resonance excitations, p + N B* , the model 
for the in-medium p propagator includes the dressing of 
its 2-pion cloud via nucleon- hole (NN^^) and delta- hole 
(AA^~^) excitations. When translated into contributions 
to the jN — > pN cross sections, the pion cloud dress- 
ing corresponds to t-channel pion exchanges (with cither 
nucleon or delta final states), as well as contact terms 
dictated by gauge invariance. The cr, /2, and Pomeron 
exchanges are not included in the in-medium p spectral 
function, where they would most notably contribute at 
large 3-momentum of the p meson. Since the outgoing p 
momenta in the CLAS kinematics are typically in a range 
around q ~ 1.5-2 GeV, it is worthwhile to estimate by 
how much the broadening in the p spectral function may 
be underestimated at these momenta. To do this, we 
neglect nuclear Fermi motion and employ the schematic 
expression relating the total absorption cross section to 
the imaginary part of the in-medium p selfenergy, 



Tp = -ImJlp/mp ^ Vrci Qn (jI°n- (1) 



Assuming a (non-resonant) pN cross section of 40 mb for 
a. q = 2 GeV p imparting on a nucleon at rest, Uiei = 0.9, 
and using qn = 0.5^0, we find Tp ~ 60 MeV. This is to 
be compared to the in-medium p-meson width obtained 
for the model employed here, amounting to ^40 MeV (for 
q = 2 GeV, qn = O.Sgg, see, e.g., Fig. 4 in our previous 
article [l^]). This estimate suggests that we may under- 
estimate the p width at high momenta by up to 20 MeV.^ 

With reasonable control over the elementary p- 
production amplitude on the nucleon, we applied the 
above-constructed model to dilepton invariant-mass 
spectra off nuclei as measured by the CLAS Collabora- 
tion. The basic assumption in this approach is that the 
dilepton production process can be "factorizcd" into a 
2-stcp process, namely, p production on a single nucleon 
followed by its in-medium propagation and final decay. 
The pertinent mass-differential cross section for a fixed 



^ However, smaller momenta contribute and the "average" density 
probed by the CLAS experiment is ~0.4po for ^®Fe targets; see 
also later discussion.. 
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FIG. 1: (Color online) Probability distributions (normalized to one) of the nuclear density at the p-meson decay points, 
dPdccs.y/dgN', left panel: using an input photon-energy spectrum as in CLAS for two differrent nuclei; right panel: using a ^-^^U 
target for different photon energies. Vertical lines indicate the bins used for the histogram. 



incoming photon energy is then given by 



dM 



(2^)2 



d^p d^q d^p' e^g^ (-IniE;'^f^^^^_(g)) 

(2^)3 (2^)4 (2^)4 (fc2)2 ^4 2 
X 5{q^ - A/2) 5(p'2 „ S\k +p-q-p')Q [kp - \p\) 
xQ{\p'\-kF) T^ik,p,q) {T^k,p,q))^ 

X {p;:M + P^M |Gj(<?)r} , (2) 

with k = {E^,k), p = {Ep,p), p' = {Ep,,p'), q = {qo,q) 
being the 4-niomenta of the incoming photon, incoming 
nucleon, outgoing nucleon (all on-shell) and outgoing p 
(i.e., dilepton), respectively. Furthermore, is the p 
photoproduction amplitude off the nucleon, G^'^ are the 
longitudinal and transverse electromagnetic correlators 
encoding the in-medium p propagator, ^ and P^'^ are the 
standard 4-dimensional projection operators. As is usual 
in many-body theory, the in-medium width (imaginary 
part of the selfenergy) of a particle is calculated in terms 
of the imaginary part of the forward scattering amplitude 
on particles of the medium, corresponding to absorption 
in the forward direction. These processes include finite- 
angle elastic scattering. Thus, Eq. ([2]) does not include 
the dilepton decay of p-mesons after being scattered elas- 
tically at finite angle. There are several reasons why this 
contribution can be expected to be small. For example, 
the nucleon on which the scattering occurs is either Pauli- 
blocked (not for forward scattering!) or ejected from the 



nucleus (i.e., corresponding to an extra particle in the 
final state). Formally, one would have to calculate such 
a contribution by inserting an extra p-N scattering am- 
plitude into Eq. together with an extra folding over 
a nuclear Fermi distribution of the struck nucleon. Thus, 
this process is formally of higher order in density. Given 
the fact that the in-medium part of the p width is basi- 
cally linear in density (e.g., cf., Fig. 11 in Ref. Q), higher 
orders due to elastic finite-angle scattering should be sup- 
pressed as well. Of course, ultimately, one should com- 
pare our predictions for total dilepton production cross 
sections to absolutely normalized data, which, unfortu- 
nately, are not available at this point. 

With a realistic photon input-spectrum consisting of 
six bins of photon energies, E^ = 1.0, 1.5, 2.0, 2.5, 3.0 
and 3.5 GeV with relative weights of 13.7, 23.5, 19.3, 
20.1, 12.6, and 10.9%, respectively, a fairly good de- 
scription of the excess mass spectra^ for deuteron and 
iron targets has been found 22|. For the full in- medium 



propagator of the p encoded in Gp the input density 
has to be specified. In Ref. [l^] this was done rather 
schematically by evaluating an average density at the 
p's final decay point, implying that the p-meson prop- 
agator quickly relaxes to the surrounding medium. The 
(average) density at the decay point was estimated as- 
suming (i) a Woods-Saxon profile for the nucleus, (ii) 
a uniform sampling of the nucleus for the production 
point (i.e., the incoming photon reacts equally likely with 
each nucleon), and (iii) a p- meson path length, collinear 
with the incoming photon momentum and estimated us- 



^ For example, Gp = {m^/gp)Dp in the vector dominance model 
(VDM). The present approach accounts for corrections to VDM 
in the baryon sector as described in Ref. |23|| . 



An "excess" mass spectrum is defined as the signal spectrum (af- 
ter combinatorial background subtraction) with the (narrow) ui 
and (f> peaks subtracted. The remaining signal is then associated 
with the p-decay contribution. 
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ing an average p momentum of (7 ~ 2 GeV (determin- 
ing its velocity), and an average in- medium p lifetime of 
Tp = i/f"i'=d = 1/200 MeV =lfm/c. The resulting aver- 
age density for ^^Fe turned out to be ^O.Sfpo: but the vari- 
ations in the normalized spectra are small when varying 
this density within ±0.1^?o (Note that, at these densities 
and 3-momenta, the width of the microscopic spectral 
function is indeed around 200 MeV, which is consistent 
with the average lifetime used for the density estimate.). 

Here we improve upon the above procedure in several 
respects. First, the width, and thus the decay time, of the 
p meson depends on its 3-momentum. We replace the av- 
erage value used before with the explicit momentum de- 
pendence given by our microscopic model for the prop- 
agator. In particular, slow p mesons decay faster than 
fast ones that consequently travel farther. This could be 
significant in view of the rather large range of photon 
energies (between 0.6 GeV and 3.8 GeV) as used by the 
CLAS Collaboration [23]. Second, the assumption that 
all p mesons of a given velocity decay after a fixed (av- 
erage) travel distance, L, is subject to corrections. Here, 
we include the exponential decay characteristics as 

dN 

— =exp(-r[eA,(f),z;]i) , (3) 

implying that even for identical kinematics and produc- 
tion point the path length varies and thus different den- 
sities at the decay point are probed. The initial creation 
points of the p meson arc still distributed according to a 
realistic density profile for each nucleus [2^1 (weighted by 
volume). Following this procedure we obtain the decay 
points and thus a distribution determining how many p 
mesons decay at a given density. The final spectrum then 
follows using this distribution in the electromagnetic cor- 
relator (rather than an average density). 

In Fig. [1] we show the resulting distributions of p-decay 
points as a function of nuclear density, for a small and 
a large nucleus integrated over the initial photon en- 
ergy spectrum as specified above (left panel), and for 
a heavy nucleus resolved into different initial photon en- 
ergies (right panel). In both distributions we observe a 
clustering of the decay points at the low and high ends of 
the available densities. In particular, the absolute max- 
imum shifts from the lowest density bin (p^ ~ 0) for 
small nuclei to the highest one {g^^ ~ Qq) in heavy nu- 
clei, and similarly when going from high to low photon 
energies. The two-hump structure of the distributions 
reflects the rather sharp transition in the nuclear density 
profiles from a dense interior to the vacuum, with a rel- 
atively thin surface layer. For the CLAS input photon 
spectrum (which ultimately determines the p's velocity 
and the time dilation of its decay time), the use of a heavy 
nucleus (^'^**U) compared to a light one (^'^C) makes a big 
difference, by increasing the interior decay fraction by a 
factor of ~3 while decreasing the outside decays by a 
factor of ~2 (see left panel of Fig. [T]). Reducing the pho- 
ton energy is particularly efficient in suppressing outside 
decays (see right panel of Fig. [T]), while the exponential 



decay characteristics induces a good fraction of interior 
decays for relatively fast p mesons. The asymmetries in 
the densities of the decay points are the main difference 
compared to the use of an average density in our earlier 
work (2^ , which is at the origin of the slight variation in 
the results shown below. 



III. SIGNATURES OF MEDIUM EFFECTS 

In Sec. Ill A wc first discuss the updated inclusive 
dilepton mass spectra in comparison to the CLAS data, 
as following from the improved evaluation of decay-point 
densities as described in the previous section. In the sub- 
sequent three subsections we then scrutinize the origin of 
the medium effects to identify observables that can im- 
prove our understanding of the medium effects and test 
consistency of spectral shapes and absolute yields. 



A. Update of Invariant-Mass Spectra and 
Comparison to CLAS Data 

In Fig. [2] we compare our earlier calculations for the 
^^Fe target with our improved calculation and with the 
CLAS data at JLab [H, H!]. It also includes the com- 
parison for the ^^C target. The average densities of 0.3gio 
for carbon and 0.4gg for iron correspond to the averages 
obtained from the density distributions in the improved 
scheme (as wc will sec in Sec. IIIIDI below, such averages 
are also estimated from the transparency ratios for total 
cross sections). One finds that the inclusion of the ex- 
plicit density distribution leads to a small downward shift 
of the p-resonance peak in the spectrum, which improves 
the description of the data compared to the calculation 
using an average density. The reason for this effect is that 
p mesons with large momentum, which predominantly 
figure into the high-mass part of the spectrum, now effec- 
tively decay at lower densities, implying reduced medium 
effects, that is, a narrower distribution (and a slightly 
smaller upward mass shift). This suppresses the high- 
mass tail of the dilepton spectrum and increases the peak 
height. Overall, the full model gives a good description of 
the data for both nuclei and provides a reasonable base- 
line for an advanced analysis of in-medium effects. We 
note that our findings arc in good agreement with the re- 
sults of the Boltzmann-Uehling-Uhlenbeck (BUU) trans- 
port simulations by the Giessen group [H, [3, [1^ based 
on the spectral function of Ref . [l2| . Also in this calcu- 
lation, a broadening of the p-meson spectral function by 
~70MeV was found without significant mass shift. 

Next, we address the influence of different model com- 
ponents on the shape and magnitude of the invariant- 
mass spectra. In Fig. [31 we display the results of the 
calculations in which (i) the resonances in the produc- 
tion amplitude are switched off (dash-dotted lines), (ii) 
the resonances in the production amplitude arc kept as in 
the vacuum (i.e., without in- medium broadening as spec- 
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FIG. 2: (Color online) Excess mass spectra for ^^C and ^''Fe (upper left and right panel, respectively) and ratios of these 
spectra to the deuterium spectrum (lower panels; in the ratios, both numerator and denominator are normalized to the 
corresponding target-mass number. A). The full calculation using explicit density distributions (solid lines) is compared to 
simplified calculations using a single average density. 



ified in Refs. [H, [131; dashed lines), and (ill) the vacuum 
propagator for the p-meson is employed (dotted lines). 
In the left panels of Fig. |31 the integrated yield for each 
scenario is normalized to the experimental data (which 
to our knowledge arc not absolutely normalized) using a 
fit.^ In the right panels of Fig. [31 only the full calcu- 
lations (sohd lines) are (x^-) normalized to the data and 
the same normalization factor is then applied to all other 
curves, which maintains the relative normalization of all 
theory curves. The following observations are made. The 
impact of the in-medium broadening of the baryon reso- 
nances on the production amplitude is essentially negli- 
gible, in both the shape (left panels) and absolute magni- 
tude (right panels) of the final dilepton spectra. Switch- 
ing off the resonance contributions in the production am- 
plitude altogether (dash-dotted lines) leads to a slight 



* Alternatively, one could normalize to the integrated strength of 
the spectrum. The differences in both procedures are negligible. 



narrowing of the spectral shape (left panels) but, more 
significantly, reduces the absolute yield of the cross sec- 
tion by more than 50% (see dash-dotted lines in the right 
panels of Fig. [3]) . This can be easily understood since the 
resonances play a significant role in the total p pliotopro- 
duction cross section on the nucleon for photon energies 
below i?^ ~ 2 GeV [l^. An important question con- 
cerns the sensitivity of the spectra to the medium effects 
in the p-meson propagator. We investigate this question 
by replacing the in-medium propagator with the vacuum 
one (resulting in the dotted lines in Fig. |3]). The spec- 
tral shape narrows, but not as much as one may have 
expected. In fact, when normalizing the theory curves to 
the experimental data, the resulting spectral shapes for 
vacuum and in-medium correlators are not much differ- 
ent for both carbon and iron targets, see left panels of 
Fig. [21 More quantitatively, for the carbon target, the 
per data point is close to 1 for both in-medium and vac- 
uum p propagator; that is, no significant distinction can 
be made. For the iron target, we actually have a prefer- 
ence of 1.08 over 1.55 in the when using the in-medium 
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FIG. 3: (Color online) Calculations for excess mass spectra off ^^C (upper panels) and ^"Fe (lower panels) targets compared to 
the CLAS data [l5l.[Ta] at JLab. The full calculations (solid lines) are decomposed by modifying different theoretical ingredients, 
by either (i) switching off the baryon-resonance contributions in the elementary production amplitude (dash-dotted lines) or 
(ii) switching off the in-medium widths of the baryon-resonance contributions in the elementary production amplitude (dashed 
lines), or (iii) replacing the in-medium electromagnetic correlator, Gp, with the vacuum one (dotted lines). In the left panels 
all spectra are (x^-) normalized to the data, while in the right panels we keep the normalization of the full model also for the 
other scenarios. 



propagator. This is in line with the analyses carried out 
in Ref. [l^ based on BUU calculations [H, [26l w ith the 
in-medium broadened spectral function of Ref. [12| , or on 
direct fits to the spectra with a Breit-Wigner ansatz for 
the spectral function [l^. We recall that the absolute 
yields of the spectra markedly increase for the vacuum 
spectral function compared to the in-medium one (which 
is more systematically quantified in Sec. IIII D]) . 

The question we want to address in the following is how 
we can better discriminate medium effects from the ex- 
perimental data. For example, for the many-body spec- 
tral function employed in our approach the broadening 
is much more pronounced at low 3-momenta (relative 
to the nuclear rest frame). The moderate broadening 
observed in the CLAS data is accounted for by the re- 
duced medium effects at large 3-momentum (following 
from the relatively soft hadronic form factors as deter- 
mined in Rcfs. [If, H^), as well as the limited fraction 



of in-medium decays even for the iron target. We inves- 
tigate several ways to better pin down medium effects. 
First, one could use heavier nuclei. This would allow 
more p mesons to decay inside the nucleus and thus aug- 
ment the medium effects (recall the left panel in Fig.[T|). 
In addition, the density dependence of the spectral func- 
tion could be studied. Second, one could concentrate 
on Icpton pairs with low 3-momcntum which should en- 
hance the decay fraction inside the nucleus because of 
the smaller distance they travel from production to de- 
cay point^. A breakdown of the data in several momen- 
tum bins would illuminate the 3-momentum dependence 



^ This poses formidable experimental challenges since in the 
present CLAS set-up only a detection of pairs with momenta 
above ^^1 GeV is feasible 
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of the spectral function, a critical property as indicated 
above. Third, as the presumably most straightforward 
option on the experimental side, one could utilize the 
absolute magnitude of the total cross section in terms of 
the so-called nuclear transparency ratio. 

The above-listed three possibilities are elaborated in 
the following three subsections IIIIBl IIII C[ and IIII Dl re- 
spectively. To facilitate direct comparisons to (existing 
or future) CLAS data, we use the same folding over the 
incoming photon-energy spectrum as quoted above p^ . 
unless stated otherwise. Our hope is that several long- 
standing issues about basic properties of the p-spcctral 
function in nuclear matter (see, e.g., the differences be- 
tween the spectral functions of Refs. [ifl,[lJ,[l3,[23|) could 
be resolved by detailed studies of accurate e+e" spectra 
in reactions off nuclei. One of our main goals is to quan- 
tify the notion of "accurate" . 



B. A Dependence of the Spectral Shape 

A systematic study of the A dependence of the abso- 
lutely normalized dilepton mass spectra under the CLAS 
conditions is compiled in the left panel of Fig. |4l One 
clearly recognizes the enhanced sensitivity to the medium 
effects with increasing A of the nucleus, mostly due to the 
simple fact that a larger nuclear volume increases the 
fraction of in- medium decays (recall also Fig. [T]). The 
"gain" in medium effects is quite appreciable when going 
from the currently available iron target to Au, Pb or U, 
which is more pronounced than the difference between C 
and Fe. It instructive to compare the ^-dependence of 
the spectra to the density dependence in infinite nuclear 
matter, as shown in the right panel of Fig. U) This com- 
parison suggests that for A ~ 200 targets the inclusive 
spectral shape (or net broadening) corresponds to rather 
moderate average densities, ~ O.S^q. In an attempt 
to better quantify the sensitivity of the spectral shape 
to the in- medium broadening, we compare in Fig. [5] the 
full calculations (solid lines) with calculations using the 
vacuum spectral function, where the latter are normal- 
ized such that the peak height coincides with the full 
calculation for each of the four nuclei. The sensitivity to 
medium effects may then be defined as the enhancement 
of the full over the vacuum curve in the low-mass tail, 
e.g., at M = 0.5 GeV. The pertinent low-mass enhance- 
ment factors are 1.35, 1.49, 1.6, and 1.73 for A = 12, 56, 
120, and 238, respectively. 



C. 3-Momentum Dependence 

Next we investigate the option of applying 3- 
momentum cuts on the outgoing lepton pair (which in 
our approximation is equivalent to the 3-momentum of 
the p relative to the nuclear rest frame). From general 
considerations (and our discussions in Sec. |ll| we expect 
that low-momentum dilepton samples should contain no- 



ticeably larger fractions of p decays in the nuclear volume. 
Figure [5] contains the results where the 3-momentum of 
the outgoing p is restricted to lie either below 0.5 GeV 
or between 2.0 and 2.5 GeV. We refer to the former as 
"low-momentum cut" , while the latter represents a typ- 
ical momentum range resulting from the photon input 
spectrum used by the CLAS Collaboration (this range 
of "intermediate" momenta accounts for roughly 25% of 
the total yield for the full calculation). The potential 
benefits of a low-momentum cut are twofold. In addition 
to the already mentioned increase of decay probability in 
the nuclear volume (which also augments the decays at 
higher densities), the medium modifications of our p spec- 
tral functions increase appreciably toward smaller mo- 
menta. The numerical results displayed in Fig.[S]corrobo- 
ratc the discrimination power of the low-momentum cut. 
The shape of the resulting mass spectrum is markedly 
different when the full calculation (dashed lines) is com- 
pared with a calculation where the vacuum propagator is 
used instead (dot-dashed lines), already for rather light 
nuclei. The main effect is a quenching of the p peak 
at its resonance position, while the (absolute values of 
the) low-mass cross section turn(s) out to be very simi- 
lar (the absorption is approximately compensated by the 
spectral broadening). Again, one may try to quantify 
the sensitivity to the medium effects by defining an ap- 
propriate ratio. The coincidence of the low-mass cross 
section suggests the peak ratio as a suitable measure, 
that is, the ratio of the resonance peak for the vacuum 
calculation to the value at the same mass for the full 
calculation. For the low (high) momentum cut we find 
ratios of 2.4 (1.2), 2.8 (1.4), 4.2 (1.5), and 4.4 (1.6) for 
A ~ 12, 56, 120, and 238, respectively, which provides a 
much more increased sensitivity over the inclusive shape 
analysis. Indeed, when comparing the cross sections for 
vacuum and in-medium spectral functions for the above 
defined "representative" 3-momentum range, the sensi- 
tivity is substantially reduced. 

D. Transparency Ratio (Yield) 

Finally, we turn to our third option to scrutinize 
medium effects, by studying the absorption of p mesons, 
rather than the spectral shape of the dilepton invariant- 
mass spectra. The width figuring into the mass spectra is 
of the same physical origin as the forward absorption of p 
mesons when traveling through the nuclear medium. Ex- 
perimentally, the nuclear absorption is usually quantified 
by the nuclear transparency ratio, 

A ■ cFyN -y px J dM {da/ dAf) vac 

which measures the suppression of the nuclear cross sec- 
tion per nucleon, relative to the cross section on the free 
nucleon (note that the latter cross section includes the 
finite decay width of the p in vacuum). This quantity 
was previously measured in nuclear photoproduction of 
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FIG. 4: (Color online) Absolutely normalized cross sections for excess mass spectra for several nuclei (left panel) from carbon 
(uppermost curve) to uranium (lowest curve) and for several fixed nuclear densities (right panel) from O.lgo (highest curve) to 
Qo (lowest curve). 



uj and (f) mesons |28l - l30l | . from which large in- medium 
widths of ~ 100-150 MeV have been extracted (see also 
Refs. m, dllll^l for theoretical evaluations). These val- 
ues are more than an order of magnitude larger than 
the respective vacuum widths of oj (F^^ =8. 5 MeV) and 
(F0=4.3MeV). This is the main reason why for narrow 
vector mesons the transparency ratio is much more sen- 
sitive to medium effects than invariant-mass spectra (in 
which the in-medium to vacuum decay fraction is small) . 
The transparency ratio is less favorable for the p meson 
where the large vacuum width suppresses the contribu- 
tion from vacuum decays (which, in principle, is of course 
a welcome feature) . In addition, the integration ranges in 
the transparency ratio, Eq. become more ambiguous 
for a broad resonance, at least in practice (i.e., in exper- 
iment). However, we already deduced from Fig. |3] that, 
also for the p meson, the change in the absolute yields 
may be more pronounced (i.e., easier to observe) than 
the medium modifications of the spectral shape, at least 
for relatively small nuclei and toward higher momentum. 
We therefore computed the nuclear transparency ratio 
for the p meson within our approach^ by integrating the 
invariant- mass spectra as computed above over M = 0.5- 
1.1 GcV. The results are displayed in Fig. [6] as a func- 
tion of nuclear mass number, A. Even for small nuclei a 
considerable suppression occurs which gradually becomes 
stronger for larger nuclei until a ratio of about 0.6 is 
reached for uranium. This is to be compared to the mea- 
surements for the and w, where the suppression reaches 
down to ^^0.4, even when normalizing to carbon. Again, 
the reason for the apparently weaker effect for the p is 
its large vacuum width of ~140McV, relative to which 



changes in the medium are at most a factor of 2-3. On 
the contrary, for the uj and (j) more than a factor of 10 
broadening was extracted and/or predicted [ill l33l - [38l |. 
which leads to much stronger effects in the transparency 
ratio. For completeness, we also investigated different 
components in our model calculation. When leaving out 
the in- medium width increase of the baryon resonances in 
the production amplitude, Ta exhibits a small enhance- 
ment of about 0.025-0.05 owing to a small increase in 
the production cross section. This effect was essentially 
invisible at the level of the invariant- mass spectra and un- 
derlines the sensitivity of the transparency ratio to even 
small variations in the calculation. The main suppression 
mechanism is, as expected, caused by the absorption of 
the propagating p meson, as illustrated by using the vac- 
uum propagator in which case Ta is close to one (the 
small residual suppression being caused by the medium 
modifications in the production cross section, essentially 
due to the in- medium width of the resonances) . 

It is instructive to evaluate the nuclear transparency 
ratio as a function of nuclear density in infinite nuclear 
matter^ (cf. Fig. [7]) . This enables yet another angle at 
the question of what typical ( "average" ) densities might 
be representative for a given nucleus. For example, we 
find that to obtain the same Ta for ^^C (^^Fe) an infinite 
density of about 0.2-0.3^)0 (0.3-0.4£io) should be used, 
which is consistent with the values obtained with the 
method used in Rcf. [23 |. 

Of further practical interest is the dependence of Ta 
on the incoming photon energy. This information is dis- 



Note that our spectra shown above were already normalized per 
nucleon |22| so that no extra factor A should be included. 



^ Note that we apply the notion of "infinite matter" only to the 
strong-interaction parts of the calculation, i.e., we still pretend 
that the medium remains transparent to the electromagnetic 
probe so that dileptons can leave without re-interactions. 
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FIG. 5: (Color online) Excess dilepton invariant-mass spectra for several nuclei. The full model (solid lines) is compared 
to the calculations using the vacuum p-propagator (scaled to the peak height of the full calculation by the indicated factors). 
Furthermore, absolutely normalized calculations are shown for in-medium and vacuum spectral functions applying the indicated 
momentum cuts to the spectra. 



played for our approach in Fig. [HI which shows that 
the suppression markedly increases with decreasing pho- 
ton energy. Lower photon energies produce p mesons of 
smaller 3-momentum which in turn increases the time 
spent in the nuclear medium and thus provides longer 
durations of suppression. In fact, for incoming photon 
energies below ~ 1.5 GeV, directed onto A ~ 200 
targets, Ta has approached the infinite-matter value for 
saturation density (0.4) within about 20%. This reit- 
erates once more the power of low-momentum cuts on 
the Icpton pair. To quantify the A dependence of the 
transparency ratio in a functional form, wc fitted our 
results in Figs. [6] and [8] with the standard power-law 
ansatz, Ta = A^ . For the individual incoming photon 
energies, the extracted values of /3 are plotted in Fig. |9l 
For the full calculation using the weighted input spec- 
trum corresponding to the JLab beam we find a value 
of /3 = —0.0983. Its magnitude is significantly smaller 
than for the available uj and 4) data, but, surprisingly, 
comparable to J/'0 suppression in fixed-target p-A colli- 
sions at high energy (see, e.g., the recent compilation in 



Ref. [1^ where the exponent a = 1 + /?; however, this 
is most likely a coincidence since the physics of p and 
Jjip nuclear suppression is presumably quite different). 
Note that the overall /3-value for the incoming photon 
spectrum would correspond to an "average" photon en- 
ergy of slightly below ^ 2 GcV (the weighted mean 
of our input distribution is E^ = 2.14 GeV). We finally 
remark that especially for light nuclei like carbon a more 
complicated fit function would be required to accurately 
reproduce the calculated values of Ta- 

IV. CONCLUSIONS 

We have conducted a theoretical analysis of invariant- 
mass spectra and total cross sections of dilcptons in nu- 
clear photoproduction for incoming photon energies in 
the few-GeV regime as recently measured by the CLAS 
Collaboration at JLab. Our approach is based on an 
effective hadronic model in which we combined in a con- 
sistent way an earlier constructed elementary production 
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FIG. 6: (Color online) Nuclear transparency ratio, Ta, as a 
function of nuclear mass number, A. The full calculation is 
compared to results when switching off the medium effects 
on the baryon resonances in the production amplitude and 
when using the vacuum p-propagator in the electromagnetic 
correlator, G„. 
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FIG. 8: (Color online) Transparency ratio as a function of the 
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FIG. 9: (Color online) Fit parameter /3 used to describe the 
dependence of the transparency ratio on the photon energy. 



FIG. 7: (Color online) Same as Fig. [6] but in infinite nuclear 
matter as a function of density. 



amplitude with an in-medium p propagator. The former 
describes well the p production off the nucleon, while the 
latter was successfully applied to dilcpton spectra mea- 
sured in high-energy heavy-ion collisions at the CERN- 
SPS. 



We first improved on our earlier calculations [22| by 
constructing a more realistic density distribution for the 
decay points of the p meson which, in particular, accounts 
for the 3-momentum dependence of the nuclear path 
length in connection with the microscopic in-medium 
width in the propagator. This leads to a slight improve- 
ment in the description of the invariant-mass spectra 
measured by the CLAS Collaboration for iron targets, 
which we also extended to carbon targets. The sensi- 
tivity to medium effects in the inclusive mass spectra is, 



however, somewhat limited in the absence of absolutely 
normalized data. Calculations for heavier target nuclei 
indicate that the situation improves when increasing the 
nuclear-mass number to A ~ 200. A more powerful lever 
arm turns out to be a cut on the outgoing lepton pair mo- 
mentum. For small momenta, say, q < 0.5 GeV, the spec- 
tral shape becomes very flat with hardly any p resonance 
peak left. The reason for this effect is twofold: on the 
one hand the kinematics make slow p mesons stay longer 
in the nuclear volume and thus enhance the in-medium 
decay fraction; but dynamics cause the medium effects 
in the spectral function to increase substantially toward 
small three-momentum. It is this dynamic effect that rec- 
onciles the large medium effect observed in heavy-ion col- 
lision (average p width of ~ 350-400 MeV [§] as extracted 
from the NA60 data on low-mass dimuons) with the 
rather moderate modifications (Tp ~ 220 MeV from the 
current CLAS data) in nuclear photoproduction. This 
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clearly calls for a low-momentum measurement in the 
latter experiments, to scrutinize the 3-momentum depen- 
dence of the medium effects as well as the prevalence of 
baryon-induced modifications in the interpretation of the 
dilepton data in ultrarelativistic heavy-ion collisions. We 
furthermore found that the nuclear dependence of abso- 
lute production cross sections can be used as a rather sen- 
sitive, complementary measure for the in-medium width. 
These are conveniently analyzed using the so-called nu- 
clear transparency ratio. While the sensitivity to an 
increase in the absolute value of the width is less pro- 
nounced for the p than for the narrow vector mesons w 
and (/), we still predict a substantial reduction of this 
ratio, which should be straightforward to measure. An 
important point here is that, for the p, the yield mea- 
surements can be effectively combined with, and tested 
against, the spectral information from the invariant-mass 
distributions (which is very challenging for the narrow 
vector mesons owing to their large decay fraction outside 
the nuclear volume). 

As a further direct extension of the present work, we 
plan to study dilepton spectra in nuclear electroproduc- 
tion (via the reaction e~ A — >■ e~ e^e~A), which would 
provide further constraints by investigating the various 
mechanisms discussed in this article under modified kine- 
matics. It would also be interesting to study tt+tt" fi- 
nal states, which are more easily accessible experimen- 
tally [i^l but suffer from extra absorption factors on the 
outgoing pious thus reducing the sensitivity to the cen- 
tral densities in nuclei. It would be very illuminating 
to perform detailed comparisons to transport theoretical 
calculations. The latter, in principle, performs a bet- 
ter treatment of the off-equilibrium aspects in the prop- 
agation process, whereas the implementation of broad 
spectral functions is more involved than in equilibrium 
approaches like the one adopted here. A continued ef- 
fort on both experimental and theoretical fronts, in both 
heavy-ion and nuclear production contexts, holds excit- 
ing prospects for further progress in the understanding 
of hadron properties in the medium. 
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